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ABSTRACT 

Jets powered by high-mass X-ray binaries must traverse the powerful wind of the companion star. 

We present the first global 3D simulations of jet-wind interaction in high-mass X-ray binaries. We 
show that the wind momentum flux intercepted by the jet can lead to significant bending of the jet 
and that jets propagating through a spherical wind will be bent to an asymptotic angle ipoo- We derive 
simple expressions for as a function of jet power and wind thrust. For known wind parameters, 
measurements of ipoo can be used to constrain the jet power. In the case of Cygnus X-1, the lack of 
jet precession as a function of orbital phase observed by the VLBA can be used to put a lower limit 
on the jet power of Ljet ^ 10^® ergss“^. We further discuss the case where the initial jet is inclined 
relative to the binary orbital axis. We also analyze the case of Cygnus X-3 and show that jet bending 
is likely negligible unless the jet is significantly less powerful or much wider than currently thought. 

Our numerical investigation is limited to isotropic stellar winds. We discuss the possible effect of wind 
clumping on jet-wind interaction, which are likely significant, but argue that our limits on jet power 
for Cygnus X-1 are likely unaffected by clumping unless the global wind mass loss rate is orders of 
magnitude below the commonly assumed range for Cyg X-I. 

Subject headings: x-ray binaries, microquasars, jets 


1. INTRODUCTION 

X-ray binaries consist of a compact object, such as a 
neutron star {e.g., Sco X-1, iHiellming et al .11 199(111 or a 
blac k hole (e.g. Cygnus X -1. lOrosz et al Mnr Cygnus 
X-3, iZdziarski et al.l 1201311 . accreting mass from a com¬ 
panion. In certain spectral states, the accretion flow near 
the compact object generat e s powerful, collim ated jets 
(jMirabel &: Rodriguez! 1199^ IGallo et al.l l20^ . These 
jets appear remarkably similar to the jets produced by 
supermassive black holes (SMBH) in active galactic nu¬ 
clei (AGN) in morphology, spectral properties, and ener¬ 
getics (when set in relation to the overall energy released 
by accretion). 

The mostly featureless non-thermal spectra of relativis¬ 
tic jets limit quantitative analysis of jet properties to 
relatively coarse estimates. However, the propagation of 
jets and their interaction with the environment offers a 
very powerful way to study jet properties that comple¬ 
ments direct studies of jets themselves. The observations 
of cavities in galaxy clusters by the Chandra X-ray Tele¬ 
scope offer an example of how jet-environment interac¬ 
tions can be used to con strain the properties of large sam¬ 
ples of AGN jets (e.g. iMcNamara fc Nnlsenl 120071 and 
references therein). 

A sub-class in the study of jet-environment interactions 
involves cases where the accreting black hole launching 
the jets is moving relative to the surrounding medium. 
In AGN, when the black hole is moving with considerable 
speed with respect to the environment, the radio emit¬ 
ting jets are swept backward by ram pressure, generat¬ 
ing a bow shock ahead of the moving black hole. AGN 
jets that show bent morphology are often called “bent 
doubles” or tailed radio sources (|Begelman et al.l [19791 : 
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iFreeland fc Wilcotsll20lH iMorsonv et al.l[2?)T^ . 

It has been suggested that a subset of low mass x- 
ray binaries (LMXBs) that move through the ISM at 
high speed due to the kick velocity the LMXB re¬ 
ceived in the supernova explosion shou ld also exhib i t bow 
shocks and trailing neck structur es (iHeinz et al.l 120081 : 
iWiersema et al.ll2009t lYoon et ^I2011h 7 very similar to 
the case of bent-double AGN sources. 

However, an aspect that makes jet propagation in X- 
ray binaries fundamentally different from AGN jets is the 
presence of the companion star. The early type compan¬ 
ion stars of high mass X-ray Binaries (HMXBs) drive 
powerful winds in the vicinity of the compact object, 
which are often the source of accretion in these objects. 
The wind mass loss rates of the OB-type donor star can 
be substantial: Mwind 10“”^ — 10“^ Mq yr“^. The stel¬ 
lar wind dominates the HMXBs’ circum-binary environ¬ 
ment compared to any winds launched by the accretion 
flow. In this work, we study the dynamics and evolu¬ 
tion of the jets affected by spherical stellar winds from 
OB-type donor stars in HMXBs. A subset of HMXB sys¬ 
tems consist of a Wol f-Rayet star in orbit with a black 
hole (e. q. Gvgnus X-3 iMartf et al.ll20(Hl: IZdziarski et al.l 
1201311 : the mass-loss rate from Wolf-Rayet stars is even 
higher than that from OB-type stars. Moreover, the 
Wolf-Rayet star in Cygnus X-3 is tidally locked to the 
4.8-h orbital period, resulting in a co nsiderable equatoria l 
enhancement of the mass-loss rate (jvan KerkwiiklfTOQ^ . 
An investigation of non-spherical and highly-flattened 
winds is beyond the scope of this study. While we briefly 
discuss the case of Cygnus X-3 under the assumption of 
a spherical wind, we plan to discuss the effects of non- 
spherical winds in a future paper. 

HMXB winds ar e not simple. For example, they 
are l ikely clumpy (lOwocki et al.l II988I : lOskinova et al.l 
I2OI2II . like winds generated by other high-mass stars. 
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The interaction of a microquasar jet with clumpy 
wi nd medium has been studied in 3D simulations 
by iPerucho fc Bosch-RamoiJ (I2012D . which we will re¬ 
fer to as P 12 hereafter, who show that clumping 
can significantly increase jet disruption. Moreover, 
even on average, they are not spherically symmetric 
due to the gravitational focusing by the compact ob¬ 
ject and Cor iolis and centrifugal effects due to or- 
bital motion (|Friend fc Castod 119821 : iMiller et al.l 120051 : 
iHadrava fc Cechurall201^ . In addition, the wind can be 
ionized by X-rays from the accretion flow, r educing or 
elimi nating line driving and stalling the wind (jGies et al.l 
l2008fl . If the X-ray flux is high enough to excite the outer 
layer of the star, a thermally driven wind may replace the 
quenched radiatively driven wind. At the current time, 
it is not clear whether the illuminated side of the wind 
would suffer from the same line-driving instability that 
generates clumps in regular massive star winds. 

In this pilot paper, we will neglect some of the more 
poorly understood complications likely present in binary 
winds and instead treat the wind as a radiatively driven 
wind that is is otropic at the surfa ce of the companion 
(following, e.g. iCastor et al.lll97a) . This will allow us 
to isolate the fundamental differences in jet propagation 
in the presence of a wind compared to jets propagating 
into uniform medium. We will discuss the effects of some 
of the likely complications in N4.7I The likely most im¬ 
portant caveat is the potential dumpiness of the wind; in 
that sense, the global simulations presented below should 
be considered a complementary approach to the detailed 
3D jet-clump simulations presented in P12. 

In §2, we present the numerical setup and the code used 
in our parameter study of jet-wind interaction. In §3 we 
discuss the results of the simulations. In §4, we compare 
the numerical results with analytic expressions derived 
for the limiting case of small deflection angles and apply 
our model to the HMXBs Cygnus X-1 and Cygnus X-3. 
In §5 we summarize our results. 

2. TECHNICAL DESCRIPTION 
2 . 1 . The FLASH Code 

Simulations were p erformed with the FLASH 3.3 
hydrodynamics code (jFrvxell et al.l l200flf l. which is a 
Message Passing Interface(MPI)-parallelized, modular, 
block-structured adaptive mesh refinement code. We 
employ the non-relativistic unsplit mesh solver, which 
solves the Riemann problem using an unsplit staggered 
mesh scheme on a three dimensional Cartesian grid 
(jLee fc DeanellW)^ . 

2 . 2 . The Wind and Jet Nozzles 


evacuation of zones adjacent to the nozzle due to the 
large velocity divergence along the jet axis. 

The stellar wind nozzle is modeled as a spherical 
boundary with inflow boundary conditions matching the 
desired wind parameters. We evolve the simulation with 
only the stellar wind present for one full orbital period 
to establish a stable, self-consistent wind profile before 
switching on the jet nozzle. 

The equation of state is assumed to be adiabatic, track¬ 
ing two separate phases of the fluid (each represented by 
a separate passive tracer fluid to distinguish wind and jet 
fluids during the computation and in post-processing). 
The wind gas is assumed to be a monatomic ideal gas 
with an adiabatic index of 7 = 5/3. The internal com¬ 
position of jets is currently unknown, however, it is rea¬ 
sonable to assume that they are strongly magnetized and 
that a sizeable fraction of their internal energy is carried 
by relativistic electrons, given the observed synchrotron 
radiation. We investigate jets composed of fluids both 
with relativistic equation of state with 7 = 4/3 and with 
a cold-gas equation of state with 7 = 5/3 and present 
results in terms of an unspecified value of 7 wherever 
possible. 

Most simulations were performed using 7 = 4/3, rep¬ 
resenting the equation of state for a gas with relativistic 
inter nal pressure, either from a fully tangled magnetic 
field (iHeinz fc Begelm^ l2000fl or a relativistic compo¬ 
nent of the gas; since our simulations are sub-relativistic, 
a relativistic equation of state implies that the inertial 
density is dominated by cold particles (e.g., protons). A 
non-relativistic equation of state represents a jet with 
internal pressure dominated by non-relativistic thermal 
plasma. As we will show, our results are only moderately 
sensitive to the actual value of 7 . 

The gravitational fields of the black hole and the com¬ 
panion are modeled as point source potentials. 

To allow for direct comparison with the analytic for¬ 
mulae we derive in fl4.2l most of the simulations pre¬ 
sented in this paper do not include wind driving by ra¬ 
diation pressure. Instead, we assumed an asymptotic 
wind at the injection at the stellar surface, i.e., a wind 
with terminal velocity Vcxi and at fixed mass flux of 
MwinH ^ 10~^ M f;i yr~^, which is typical for OB stars 
(iPuls et al.ll2008ll . 

Wind driving was incorporated in a sub-set of the sim¬ 
ulations presented in this paper to test the sensitivity of 
our results again st th e assumption of an asymptotic ra¬ 
dial wind. (See 114.31) . For the radiatively driven wind, 
the initial velocity field follows the so-called /3-law 

i^(r) = noo(l - ro/r)^, (1) 


Both the jet and the wind injection are modeled as 
inflow-boundary conditions on an interior portion of the 
grid (we will refer to the regions excluded from the hydro- 
dynamic integration and instea d treated as an in terior 
boundary as “nozzles” following iHeinz et al.ll2006ll . 

The jet nozzle has a cylindrical shape with inflow 
boundary conditions at the surface, injecting a bipolar 
outflow with a prescribed energy, mass, and momentum 
flux to match the parameters we choose for the jet. For 
reasons of numerical stability, we inject a slow lateral 
outflow from the side walls of the cylinder with negli¬ 
gible mass and energy flux in order to avoid complete 


where v^c = 2, 500 kms“^ is the wind terminal speed 
[chose n to match the w ind parameters of typical OB-type 
stars (jPnls et al.l 120081 )] and rg = i?*x, where i?* is the 
radius of the star and x is applied to avoid zero velocity 
and infinite density on the surface of the star. We set 
the stellar radius i?* = 1.4 x 10^^ cm. The value of x 


can be expressed as x 


1 - (u*/?;oo)^^^ 


= 0.99, where 


the wind surface velocity, v^,, is of the order of 10 ®cms“^, 
which is the sound speed with Tgs ~ 30, 000 K. However, 
due to the steep variation of density and pressure around 
the surface, there is a limit in performing a numerical 
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calculation with such a high value of x. Alternatively, we 
set the value of a: = 0.95, which was chosen to be high 
enough to maintain the initial wind profile from a typical 
OB-type star by iterative ID simulations. The mass flux 
of the wind was fixed at the surface of the star (where 
density and velocity of the injected wind determine M 
uniquely) and li ne driving was modeled using the Sobolev 
approximation (iCastod Il974f) , such that the total line 

/ , , \acAK 

acceleration results in oc ( p 'TF j ; where acAK 

is the parameter of the CAK model ([Castor et al.lll97^ . 
typically depending on the effective temperature of the 
star. We choose a value of acAK = 0.64 in our model 
corresponding to a typical OB-type star. 

Our simulations are adiabatic and scale-free. In phys¬ 
ical units chosen to approximately match the wind 
and binary properties of Cygnus X-1 and allow simu¬ 
lations to be completed within the available computa¬ 
tional resources, the jet velocity was set to be ujet = 
3 X 10®cms“^, with an initial internal Mach number of 
A4jet,o = 30 at the base of the jet (the “nozzle”). To 
explore the dependence on Mach number, we ran a sim¬ 
ulation at Aljet.o = 10 and otherwise identical param¬ 
eters compared to our fiducial run with a jet power of 
Ljet = 10^® ergss“^. Note that the Mach number Aljet 
varies along the jet given the adiabatic behavior of the 
fluid. 

In order to resolve the hydrodynamics at the injection 
scale with at least 10 cells across the nozzle, we forced 
the jet nozzle to be at maximum refinement, resulting 
in an effective resolution of 4.7 x 10® cm = 1.6 x 10“^ a, 
compared to an orbital separation of a ~ 3 x 10^^ cm for 
Cygnus X-1 (|Gies &: Boltonl[l9^ . The radius of the jet 
nozzle is about 2.5 x 10^® cm and the full box size of the 
simulation is about 4 x 10^^ cm, centered on the center 
of mass of the binary. 

We varied the jet power to span the range Ljet ~ 
10 ®®, 10 ®®, 10 ®^ ergs s“®, comparable to the range of un¬ 
certainty in the jet power of Cygnus X-1, 9 x 1 0®® — 
10®^ergss“® (jCailo et al.l 120(151 : iRussell et al.ll2007[ l. 

For the bulk of our simulations, the jet was injected in 
a direction perpendicular to the orbital plane, le., along 
the z-axis of our grid. We also investigated off-axis jets 
with angles of 30°, 60°, 75° relative to the orbital axis of 
the system, inclined towards the binary companion (in¬ 
clination angles perpendicular to the orbital separation 
vector would not result in any change in the simulation, 
given that simulations only cover a small fraction of the 
binary period once the jet is switched on). The detailed 
model parameters of our different runs are described in 
table [TJ 

Note also that in cases where the jet is oriented per¬ 
pendicular to the binary separation vector a, the sim¬ 
ulations have mirror-symmetry about a; in these cases, 
in order to reduce the need for computational resources, 
we only simulate the upper hemisphere at full resolution, 
but include both hemispheres to avoid spurious bound¬ 
ary effects near the orbital plane. Results in those cases 
are quoted for the high-resolution half of the simulation. 

2.3. Orbital Motion 

The binary parameters were set loosely approximate 
the parameters for of Cygnus X-1. For simplicity, we set 


the mass of the black hole and the star to be 10 and 20 
M 0 , respectively, and the separation between them is set 
to be 3 X 10^^ cm, which gives an orbital period of 5.8 
days for Cygnus X-1, compared to the observed val ue of 
5.6 days (iBrocksoop et al.lll999HPoolev et al.lll999[l 

As we show below, the jet propagation time across one 
binary separation (1.7 minutes) and the time required 
for a quasi-stationary bent jet solution to be established 
(approximately 10 hours) is much shorter than the or¬ 
bital time. For numerical simplicity and to allow direct 
comparison with the analytic formulae presented in 94.21 
we neglected orbital rotation in most of our simulations, 
keeping the two nozzles stationary in our cartesian grid. 
As a first order approximation, this is justified because 
the orbital velocity is only of order 20 % of the wind ve¬ 
locity, and thus orbital effects on the wind ram pressure 
introduce corrections of the order of only 5%. 

In order to verify that the effects of orbital motion on 
the gross dynamics of jet propagation are small, we ran 
a sub-set of the simulations including orbital rotation, 
presented in 94.21 In this case, the nozzles (star and jet) 
move along their orbital trajectories in the x-y plane (be., 
we simulated the orbit in a fixed, non-rotating frame, 
which eliminates the need to introduce terms for Coriolis 
and centrifugal forces into the solver). The coordinate 
origin was set to be at the center of mass. 

In the rotating case, we assumed that the star is co¬ 
rotating with the orbit, such that the outflow velocity at 
the stellar surface is given by 

"Cwind.rot (^) — CJ X X nwind,* , (2) 

where uj is the orbital angular velocity, Fwind,* the wind 
velocity at the stellar surface calculated from eq. CD, and 
i?* the position of the star. 

We ran the simulations for one orbital period before 
switching on the jet in order to allow the flow to estab¬ 
lish a converged velocity and density profile. After the 
jet launches, the simulations were carried out for several 
hours in real time units, long enough to establish the bow 
shock and the jet in a quasi-steady state (see 92.41) . 

Finally, we performed a comparison simulation of a jet 
propagating into a uniform wind with parameters match¬ 
ing those of our hducial run at the position of the com¬ 
pact object (referred to below as UniWind_E36). 

2.4. Measurement of the Jet Thickness and Propagation 

Direction 

A key variable determining the strength of the jet-wind 
interaction is the thickness h of the jet as seen by the 
wind (be., the size of the jet perpendicular to both jet 
and wind velocities). We measured h as follows. 

In post-processing, we identified matter inside a com¬ 
putational cell as jet material if the value of C = 
{vz/vjet)J exceeded a fixed threshold, where {vz/vjet) 
is the flow velocity normalized by the initial jet speed, 
which is 3 X 10®cms“® in our standard parametrization, 
and J is the fractional density of jet tracer fluid injected 
at the nozzle. We found that a choice of C = 0.1 success¬ 
fully identified the jet material in all cases (see Figure[7]). 
The thickness of the jet was measured in the direction 
perpendicular to the separation vector and the jet axis, 
since it is the dimension of the jet in that direction that 
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Table 1 

Model Parameters 


Model 

Ljetfergss l) 

■AdjetjO 

inclination (degree) 

hi^ (cms-l) 

SphWind_E35 

1035 

30 

Oa 

3 X IQio 

SphWind_E36 

1036 

30 

0 

8 X lO^o 

SphWind_E36_M10 

1036 

10 

0 

1 X lO^i 

SphWind_E36_rot 

1036 

30 

0 

8 X 103° 

SphWind_E36_acc 

1036 

30 

0 

8 X 103° 

SphWind_E36_30deg 

1036 

30 

30 

8 X IQio 

SphWind_E36_60deg 

1036 

30 

60 

8 X IQio 

SphWind_E36_75deg 

1036 

30 

75 

8 X lO^o 

SphWind.E37 

1037 

30 

0 

1.5 X 10“ 

SphWind E37 gaml66 

1037 

30 

0 

1.5 X 10“ 

UniWind_E36 

103B 

30 

0 

8 X lO^u 


0° indicates that the direction of the jet is perpendicular to the line between the star and 
the Black hole. 

The jet thickness at the re-collimation shock, hi, is measured naively by checking the 
variation of the jet thickness along the jet from simulation results. 


determines the amount of wind momentum flux inter¬ 
cepted by the jet. 

When the jet first turns on, it propagates in its initial 
direction, following the standard evolution of jet propa¬ 
gation until the expansion of the cocoon becomes slower 
than the wind velocity. As the head of the jet propagates, 
the accumulated perpendicular momentum flux begins to 
bend it away from its initial propagation direction. We 
then traced the propagation direction of the jet fluid to 
determine the jet trajectory and bending angle. 

3. RESULTS 

3.1. Jet Bending in Spherieal Winds 

Our simulations confirm the general expectation that 
a powerful wind from a companion star can affect the 
propagation of the jet, and that the ultimate trajectory 
of the jet depends on the relative momentum flux in the 
wind and the jet, as well as the geometry of the system. 

We briefly describe the morphology of the jet-wind in¬ 
teraction and compare it to simulations of jet bending ob- 
served in the inte raction of a jet with a uniform medium 
(|Yoon et al.ll2011h . Figured] shows snapshots of our fidu¬ 
cial run at a jet power of Ljet = 10^® ergss“^. 

Upon injection into the grid, the jet fluid is generally 
over-pressured compared to the external pressure and the 
ram pressure of the stellar wind (the pressure is fixed by 
our choice of the jet power, the jet velocity, the Mach 
number, and the cross sectional area of the jet nozzle). 

We chose this setup to allow the jet to establish a self- 
consistent, stable structure by letting the jet reach pres¬ 
sure equilibrium with the bow shock. During this phase, 
the jet does not experience any bending, given that it is 
strongly over-pressured with respect to the ram pressure 
in the stellar wind. We show a typical setup in Figured! 
which displays a density slice through the simulation. 
Generally, the lack of external confinement leads to free 
lateral expansion of the jet with a half-opening angle of 
order oq ^ l/Aljet,o- The wind is gravitationally focused 
by the black hole in the down stream region, generating 
the density enhancement along the equatorial plane vis¬ 
ible to the left of the black hole in Figure dl Here and 
in the following, we place the x-axis along the orbital 
separation vector at t=0 and the z-axis along the orbital 
angular momentum vector. 

The free expansion of the jet proceeds until its pressure 


reaches the pressure behind the bow shock of the wind, at 
which point the jet goes through a re-collimation shock 
and reaches a stable equilibrium thickness hi. Thus, hi 
is not directly set as a simulation parameter, but instead 
determined by the initial Mach number of the jet. 

However, since neither the re-collimation region nor 
the initial Mach number of the jet are observable, we will 
carry out most of the analysis in this paper using the jet 
thickness h beyond the re-collimation shock, relating it 
where possible to observables like the large scale opening 
angle of the jet Oobs- We will describe the details of the 
re-collimation region, which is very close to the black 
hole compared to the size of the simulation box, in 1)3.21 

The approaching stellar wind material goes through a 
stationary bow shock as it is forced to propagate around 
the jet. Beyond the re-collimation shock, the transverse 
pressure gradient imparted on the jet by the lateral mo¬ 
mentum flux of the stellar wind then gradually bends 
the jet fluid away from the companion star, while the jet 
thickness h is set by lateral pressure balance with the 
wind bow shock pressure. 

The effects of the radially declining wind density and 
the changing velocity of the wind as a function of distance 
along the jet imprint a qualitatively different asymptotic 
behavior of the jet compared to interaction with a uni¬ 
form medium. Because the density declines roughly as 
r“^, where r is the distance from the center of the star, 
the effect of bending declines with distance, and most of 
the bending occurs within about a binary orbital sepa¬ 
ration from the jet nozzle. 

More importantly, jet bending is caused only by trans¬ 
verse momentum flux, which depends on the angle = 


sm 


-1 


^jet ^ ^wind 


between the local jet velocity and 


wind velocity through sin^ {§). Because d decreases with 
increasing r, the amount of transverse momentum flux 
also decreases with r. Asymptotically, any initially large 
i? will tend to zero (even in the absence of any jet bend¬ 
ing) and the amount of lateral momentum flux across the 
jet decreases strongly with distance. At infinity, jet and 
wind travel parallel to each other at some asymptotic an¬ 
gle tpoD relative to the initial jet direction. In contrast, in 
a uniform wind, d does not decrease with distance from 
the nozzle, and both jets must eventually be bent such 
that the jet plasma asymptotically flows parallel to the 
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wind direction. 

We ran simulations with jet powers of Ljet = 
10^^, 10^®, 10^^ergss“^, and the results are shown in 
Figure [5] The dashed lines indicate the converged 
asymptotic lines towards which the jet is bent by the 
wind, showing that the ipoo is a strong function of jet 
power. In the case of Ljet = 10^® ergs s“^, the ram pres¬ 
sure by the stellar wind is sufficiently strong to bend the 
jet by almost 90°, similar to the case of the UniWind_E36 
uniform wind model. On the other hand, for the highest 
jet power case in our simulation, Ljet = 10^^ ergs s“^, the 
bending angle is small. We will discuss the relationship 
between the jet kinetic power and the inclination angle 
in 114.21 

We generally find that the jet is bent towards its 
asymptotic bending angle within a time scale of Tbend 
O.ObTorbit, where Torbit is the orbital period; this is 
roughly the time for the geometry of the inner jet to 
reach a steady state. We measured the propagation di¬ 
rection and asymptotic bending angle of the jet after it 
settled into its steady state. Note that the short bending 
time scale Tbend << Torbit implies that the jet reacts in¬ 
stantaneously to changes in binary orbit, which further 
implies that the jets must be precessing on the orbital 
period of the system if bent by jet-wind interaction. 

Figure |3] explores the dependence of our simulations 
on the initial internal Mach number of the jets, with 
A4jet,o set to one third of our fiducial value. With oth¬ 
erwise identical parameters, a smaller Mach number im¬ 
plies larger thermal energy relative to the total energy of 
jet. The increased thermal pressure leads to a larger 
initial opening angle of the jet, which in turn results 
in an increase in transverse momentum transfer and jet 
bending. The increased surface area and decreased Mach 
number also increase the incidence of Kelvin-Helmholtz 
instability and earlier onset of jet disruption. 


3.2. The Re-collimation Shock 

In order to investigate the physics of jet re-collimation 
by the wind, we carried out a set of test simulations 
with significantly increased resolution, restricted to a 
shorter duration. We performed the tests with two jet 
Mach numbers, A4jet,o = 10,30. Snapshots of the re- 
collimation region are shown in Figure SI We measured 
the jet thickness along the y-axis (y-z slice) for the anal¬ 
ysis below because the effective cross section of the wind 
momentum flux captured by the jet depends only on the 
width of the jet in y-direction. Jet bending is facilitated 
by the pressure gradient in the x-z plane, where only the 
leading edge of the jet is subject to the increased pressure 
behind the bow shock. 

In our simulations, the jet is initially freely expanding. 
Acceleration of the lateral expansion will become ineffi¬ 
cient once the lateral motion itself becomes supersonic. 
This sets the characteristic semi-opening angle ag of such 
a supersonic “fan” simply as 

1 Cs,0 nPo 1 

=- = \ -■ w 

A4jet,0 'Cjet V Po 'Cjet 


From the numerical experiment, the initial half¬ 
opening angle oq of the jet is roughly 


ao 


3 

A4jet,0 


( 4 ) 


slightly larger than the simplistic estimate oq 

1/A4jet,0- 

Once the conical expansion has been established, the 
lateral ram pressure Pjet,ram,_L of the jet 


P 


jet,ram,_L — Pjet Sin Ckg "^jet — PO 



Oo ffet 


( 5 ) 


is always larger than the internal pressure (since the lat¬ 
eral expansion is supersonic). 

In terms of the kinetic jet power 


Ljet,kin — /4 — Oq 

the lateral ram pressure is 


P 


jet,ram,_L 


jet, kin 

TTz'^Vjet 


( 6 ) 

( 7 ) 


independent of sinoo and AJjet.o- 

Lateral expansion will proceed until Ljet,ram, ± drops 
below the pressure in the wind bow shock that forms 
around the jet, L’wind,ram- At this point, a re-collimation 
shock must form in the jet and bring the internal pressure 
of the jet into equilibrium with the bow shock. We will 
denote the location of the re-collimation shock along the 
jet as zi- For parameters considered in this paper, zi is 
always much smaller than the binary separation a, so we 
will assume that ram pressure of the wind is constant for 
the discussion of zi and hi, so the wind ram pressure is 
given by its value in the equatorial plane at the location 
of the black hole. 

In terms of the mass loss rate of the winc0, Mwind = 
dTrr^riwindPwind = 47ra^UwindPwind,o, the wind ram pres¬ 
sure at the jet nozzle is then 


Lwind,ram,0 — Pwind.O^wind — 


Mwind'Cwind 

dTTa^ 


( 8 ) 


where, Pwind,o is the wind density at the footpoint of the 
jet and Uwind is the velocity of the stellar wind, assumed 
to be constant. 

The location of the re-collimation shock zi is given by 

equating Lje^ _ 1 _ = Lwind,ram,0- 


4L 


Zi = a. 


'jet 


MwindTwindTjet 


( 9 ) 


again independent of AJjet.o and Oq- 
This simple picture is qualitatively confirmed by the 
high-resolution simulations of the re-collimation shock: 
Figure [5] shows the re-collimation shock at z ~ 7x lO^^cm 
regardless of the jet Mach number, consistent within 
an error range of 20% with the analytic solution zi ~ 
9 X 10^^ cm from eq. m for the simulation parame¬ 
ters (Ljet = 10^Vgss“\ M = 10“®MQyr“\ Uwind = 
2.5 X 10® cm, 'Cjet = 0.1 c, a = 3 X lO^^cm). 


In Figure!^ we plot the measured jet thickness h as a 
function of height z. As expected, the jet with the ini¬ 
tially higher Mach number has a narrower opening angle. 


^ Given typical properties of OB-type stars, the mass loss rate is 
of order Mqb 10“^ 10“® Mq yr“^ with a terminal velocity of 

'I’wind ~ 2000 — 3000 km s“^ IlCastor et al.11137^ : IPuls et al.1120081'1 . 
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Figure 1. Time sequence of density maps for our fiducial simulation SphWind_E36. The black circle indicates the surface of the companion 
star. The enhanced density in the down stream of equatorial plane to the left of the black hole is due to the gravitationally focused wind. 
The bow shock structure along the jet reaches steady state approximately in 12 hours after it launches. The magenta are a in the right-most 
image indicates the jet materials (marked only in the lower half of the image), identified by a certain threshold (see Ej. The cyan dashed 
lines in the right-most panel indicate asymptotic lines along which the jets converge, showing that the jet is bent by approximately 30° 
from the initial direction. In the down-stream region, the bow-shocked wind passes around the jet and re-collimates in an expansion fan 
and a (weak) re-collimation shock, as expected for super-sonic flow around an object, leaving the post-shock region filled with wind gas, 
visible in the right-most two panels. 



-ixio’^ 0 -1x10'^ 0 

X [cm] X [cm] 


Figure 2. Density maps for the case of SphWind_E35 (left panel) 
and SphWind_E37 (right panel) after the steady state bow shock 
structure has been established. While the jet for SphWind_E35 is 
disrupted within a short distance from its injection, the jet for Sph- 
Wind_37 is steadily maintained (marked in magenta color). The 
cyan dashed lines indicate that the jet bending angles are 65° and 
8° for SphWind_E35 and SphWind_E37, respectively. The vertical 
black thin trajectory in right panel is the low density area gener¬ 
ated by the shear layer between the jet and the bow shock. 

At the re-collimation shock, the jet will have a thick¬ 
ness hi and will be in pressure equilibrium. Beyond ^i, 
the jet will thus adjust its thickness h to maintain pres¬ 
sure equilibrium with the wind. We will discuss the prop¬ 
agation of the jet in this phase, and the interaction with 
the wind that occurs beyond zi, in the next section. 

4. DISCUSSION 

As described above, the development of an asymptotic 
bending angle '0oo is expected from simple considerations 
of the geometry of the interaction between jet and wind. 

A full analytic description of the detailed dynamical 
evolution of the jet (including the onset of dynamical 
instabilities such as Kelvin-Helmholtz instability) is be¬ 



-1x10'^ 0 

X [cm] 


Figure 3. Density map in the case of lower Mach number, A^jet,o 
= 10 (SphWind_E36_M10). 

yond the scope of this paper. Therefore, for the follow¬ 
ing analysis, we make the assumption that the bending 
angle tp of the jet is small (implying relatively weak in¬ 
teraction) . 

Because non-linear effects (like dynamical instabilities) 
tend to increase the cross section of the jet, they will 
tend to increase the bending angle due to the large net 
transverse momentum intercepted by the jet. In that 
sense, the relations derived below will be lower limits 
on the actual bending angle. This is borne out by our 
simulations (See Figure [9] below). 

We will further assume that the wind is asymptotic, 
i.e., has reached constant velocity before interacting with 
the jet and thus follows a simple density profile. We 
will also neglect effects of orbital motion in the analytic 
approximations below (justified by the fact that they can 
be expected to be about an order of magnitude smaller 
than the dominant effects, as argued above). 
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Figure 4. Density map for re-collimating jets in A^jet,o = 10 
(upper panels), and A^jet.o = 30 (lower panels). The jet is re¬ 
collimated in the y-z plane while bending occurs in the x-z plane. 

4.1. The Evolution of the Jet Thickness Beyond the 
Re-Collimation Shock 


We will assume that the jet is in pressure equilibrium 
with the ram pressure of the wind, 


- jet,ram,_L 


Pwind(^) (^'windCOS^ 
Pwind,0^wind ^ 


Pwind,0^wind 



2 


( 10 ) 


where we have used cos 6 = aj^Ja? + 2 :^, where Q is the 
angle between the orbital separation vector a connecting 
the star and the black hole, and the vector r from the star 
to a given position along the jet. The jet has an initial 
jet thickness hi, as discussed above, which we take as a 
parameter in the following. 

The geometry of the jet and the re-collimation shock 
and the definitions of the relevant angles and coordinate 
axes are sketched in Figure El 

Beyond the re-collimation shock, the jet thickness h 
follows from pressure equilibrium between the jet and 
the bow shock: 


7jet,rani,± — Pwind.O^wind 1 2 1 2 

(X z 


— -^et — -^eq,l 


Hz) 

hi 


-27 


( 11 ) 


where Peq,i and hi are the pressure and the jet thickness 
at the re-collimation shock. This sets the jet thickness 
h: 

Figure [7] shows the measured jet thickness h as a func¬ 
tion of 2 ; for simulation SphWind_E37, compared to the 
value calculated fro m eq . (HU), for different choices of the 
jet threshold (see 112.4I) . The figures show good agree¬ 
ment between the model and the simulation. 


Generally, hi will not be measurable. However, we can 
us eq. m to relate a measured jet thickness (or an upper 
limit) at large z to the jet thickness at any other z, given 
v alues for y and m _ 

[Stirling et al.l (|200H) reported that the VLBA jet of 
Cygnus X-1 has a semi-opening angle of avLBA = 
h/2zvLBA ^ 2°, where zvlba is the scale length of the 
extended jet on which the opening angle is measured; the 
orbital separation and the length of the jet are 0.1 mas 
and 15 mas, respectively. By using eq. (fT2)l with 7 = 4/3, 
the jet thickness at the re-collimation shock can be esti¬ 
mated as hi < 5.7 X 10“^ a. For a value of 7 = 5/3, we 
find a value of hi < 1.3 x 10“^ a. We will further discuss 
the constraints on the initial half-opening angle in the 
case of the Cygnus X-1 jet in M.4I 


4.2. Jet Bending and the Asymptotic Bending Angle 

With an expression for the jet thickness h from eq. (IT2)l . 
we can now discuss the amount of bending experienced 
by the jet. Technically, jet bending occurs because a 
transverse pressure gradient exists behind the bow shock 
that the wind drives around the jet, such that the ex¬ 
ternal pressure at the leading edge is Pi ^ Pbow and 
the pressure on the trailing edge of the jet is P 2 <C Pi- 
Thus, a transverse pressure gradient exists inside the jet 
as well, acting to accelerate/bend the jet fluid away from 
the star. 

In deriving an estimate for ‘l|)oo^ we will make the sim¬ 
plifying assumption that the bending angle is small, i.e., 
that the accumulated transverse momentum flux is small 
compared to the lateral momentum flux. The reason for 
this assumption is that the jet will be dynamically dis¬ 
rupted if the bending angle is large, as the jet-boundary 
interaction must be significant in this case. Our simula¬ 
tions bear out the validity of this assumption at least in 
the hydrodynamic case studied here. Similarly, we ne¬ 
glect the momentum transfer from the wind to the jet in 
the longitudinal direction, which would lead to accelera¬ 
tion or deceleration by a small amount. 

Based on the properties of observed astrophysical jets 
(the presence of shocks, the inferred large kinetic power 
compared to the minimum internal energy based on 
the observed synchrotron intensity), we further assume 
the jet to be supersonic (large internal Mach number) 
and ballistic (no further acceleration beyond the nozzle). 
This reflects the setup of our simulations. 

Under the assumption of constant longitudinal jet ve¬ 
locity and -Mfet 0 ^ longitudinal jet momentum 

per unit jet length is conserved and given by 


^m,jet — dylj_PjetUjet 

“ ’’jet.O Pjet,0 ^^jet 
Tjet.kin 


V 


jet 


(13) 


where dA± is the area element perpendicular to the ini¬ 
tial jet direction, rjet.o and Pjet,o are the radius and the 
density of the jet at the nozzle. 

The transverse momentum per unit jet length accumu- 
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Figure 5. Left panel: Identified jets from simulation for the case of A^jet,o = 10? 30. Right panel: Measured jet thickness along jet. The 
dotted lines indicate the location of re-collimation shock. 

lated by the jet can be derived as a function of .s: 



Figure 6. Schematic figure, a is a jet semi-opening angle, a is a 
separation, 0 is an inclination angle from the orbital plane, and hi 
is the jet thickness at the collimation shock. 



Figure 7. The thickness of the jet in SphWind_E37 m odel . The 
solid black line represents the analytic solution from eq. m with 
the parameters appropriate for the model of SphWind_E37, and 
dashed and dot-dashed lines indicate the numerical results for dif¬ 
ferent choices of the threshold used to determine whether a com¬ 
putational cell belongs to the jet. 


A$ 


ft(z) 


m.wind — 


dt 


-f-l 

Jo I'jet J- 


fh/'l 

-hjl 

hjl 


dy 


rx{v)l2 


dxS/ x^ho 


l-X(y)/2 


dyPho 


0 ^jet J-h/2 

= [ —h[z')Poo^ 

Jo "^jet 

- r ^ 2 . 

— / Pwind,0^wind 

Jo "^jet 

^iPwind.O^wind cf \ 


a? + 


2-1/7 


where 


lljet 

^iMwindllwind 

47 r a i;jet 

z / a 




= f 

Jo 


dy 


1 + r 


2-1/7 


(14) 


(15) 


which can be expressed as a combination of hypergeomet¬ 
ric functions, but is most easily evaluated numerically. 

In the first order (small bending angle) approximation 
we are making here, the ratio of transverse to longitu¬ 
dinal momentum is equal to the bending angle (le., the 
angle between the local and the initial velocity vector or 
tangent vector of the jet as a function of z): 


/ / \ l^_L Adiwind 

V{Z,J) = - = —T - 

V\\ ®m,jet 

hlwindllwindlljet h 


Air a lvjet,kin 


-f{z,l) 


(16) 


The asymptotic value for z — > oo can then be evalu¬ 
ated in terms of elementary Gamma functions by taking 
the appropriate limit of f{z,j): 


I 1 . / / \ MwindliwindUiet hi , , 

i/ioo= lim i/>(z, 7 ) =----- hm /(z, 7 ) 


z^co ' 47raLjet,kin z-^co- 

Mwindllwindlljet ^1 ^r(3/2 - I/ 7 ) 

47raLjet,kin 2 r(2-l/7) 


(17) 
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where 


/(t) = lim 

z—^oo 


v^r(3/2 - 1/7) 
2 r(2-l/7) 


(18) 


with /(4/3) = 1.2 and /(5/3) = 1.07 for a relativistic 
and non-relativistic monatomic gas, respectively. Hence, 
for otherwise identical parameters, the jet should be bent 
about 10% less in the case of 7 = 5/3 compared to the 
case of 7 = 4/3, showing that the effect of adiabatic index 
on ipoo is moderate. Figure |8] shows the dependence of 
the f{j) on the adiabatic index from 7 = 4/3 to 7 = 5/3. 

While hi cannot be measured observationally, we can 
express eq. (HU in terms of the observable jet width 
on VLBA scales, /lobs measured at distance Zobs ^ 
or alternatively, the observed opening angle aobs = 
hohs/ (2Zobs) ■ 


V^oo - 


hohs 


^obs 


1/7 


Mwindl'wind'yjet 


47rL 


jet,kin 


fh) 


' ^obs 


/ ■2^obs \ ^ ^wind^wind^jet r/ \ 


27tL 


'jet,kin 


(19) 


where Zobs = Wlba/ sin(0LOs) and ctobs are assumed 
to be corrected for foreshortening given the line-of-sight 
inclination angle 0 los of the jet. 

We carried out one test simulation, Sph- 
Wind_E37_gaml66, with 7 = 5/3 for the jet fluid, 
and the resulting jet appears more straight, consistent 
with analytic expression (see Figure [S]). Note that the 
ratio is independent of the orbital separation, a, because 
of the relationship hi = 2zisinao, where zi can be 
calculated from eq. (En¬ 


in the small bending angle regime, our analytic ex¬ 
pressions are consistent with the simulation results, as 
shown in Figure |9l In the figure, the solid curves were 
constructed by integrating the jet trajectory along z, 
given the analytic expression for the transverse velocity 
■L'jet,_L ~ i’Vjet from eq. (I16|) in the small angle approxima¬ 
tion [i.e., d(x — Xh)ldz = tanip where xq is the location 
of the black hole]. For small bending angles, the figure 
shows excellent agreement between the model and the 
simulations. 

As expected, for the case of stronger jet power, Ljet = 
10^^ergss“^, the jet is only moderately bent from the 
initial direction, while lower power jets are more strongly 
affected by the winds, resulting in a higher degree of 
deflection. In the figure, the solid lines represent analytic 
trajectories calculated by eqs. (IT^ - (fT4ll and dashed lines 
indicate the asymptotic direction estimated in eq. ( 113 - 

In the case of our fiducial simulation (jet power Ljet = 
I036 ei-gss“i), the jet becomes dynamically unstable 
around z = 2.5 x 10^^ cm, which leads to significant 
broadening of the jet and enhanced bending, and, as a re¬ 
sult, the jet begins to deviate from the analytic estimate. 
Further quantitative analysis of the turbulent structures 
and their effect to the evolution of the jet is beyond this 
work. It is, however, clear from the simulations that jet- 
instability will only increase A$m,-wind/4*m,jet and '0oo. 

Thus, our small jet bending angle approximation can 
be considered a robust lower limit of the actual jet bend¬ 
ing angle. Figure [TOj shows the fractional deviation of the 
numerical result from the analytic estimate as a function 
of jet bending angle. Our analytic formula is very ac¬ 
curate for bending angles smaller than 20°. For larger 
bending angles, the approximation breaks down. 



Figure 8. The asymptotic jet bending angle as a function of adi¬ 
abatic index, 7 . The angle is normalized by the one for the case of 

7 = 4 / 3 . 

While these expressions are strictly non-relativistic, it 
is straight forward to show that in the ultra-relativistic 
case, the estimate for the bending angle "0 is increased by 
a factor of over the non-relativistic case (where the 
velocity is simply set to ?; = c). Thus, the lower limits 
we derive below on Ljet from observational upper limit 
on becomes stronger in the relativistic case. 


4.3. The Effects of Orbital Motion and Wind 
Acceleration 

We carried out two simulations including the effects of 
orbital motion(SphWind_E36_rot) and radiatively driven 
winds (SphWind_E36_acc) in order to evaluate the im¬ 
portance of both effects by comparing to our standard 
model. 

Eigure IHj shows that, for small bending angles, the 
effects are small and the results with and without or¬ 
bital motion and wind acceleration are consistent with 
each other. This is not surprising, since the centrifugal 
force from the orbital motion is about 3 orders of magni¬ 
tude less than radiative force, and the Coriolis force acts 
purely in the transverse direction. Therefore, our setup 
with fixed black hole and star positions is sufficiently ac¬ 
curate in the context of this analysis. 

We reach a similar conclusion about the effect of ra¬ 
diative acceleration. The model with a radiatively driven 
wind has a negligibly small difference in the momentum 
flux of the stellar wind at the binary separation compared 
to our standard model. It implies that our assumption 
that the wind reaches terminal velocity before encoun¬ 
tering the jet is valid. 

Eor large bending angles, where dynamical instabilities 
lead to rapid jet disruption and increased bending, the 
deviation between individual simulations is noticeable, as 
expected given the time variability of the jet trajectory 
on those scales. 
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4.4. Jet Bending as a Diagnostic of Jet Power: The 
case of Cygnus X-1 

In the limit that Afi>m,wind > ‘hmjet, bending in the 
simulation is so strong that the jet is dynamically dis¬ 
rupted, rather than simply bent. The interaction dis¬ 
perses the jet into a broad, no longer collimated flow at 
much lower velocity than the jet velocity. We would not 
expect such a strongly bent jet to survive as an observ¬ 
able radio jet. This suggests a simple diagnostic: If a 
stable compact jet is observed intact in an HMXB, one 
can conclude that the bending angle should be moderate. 

For example, the compact VLB A jet of Cyg X-1 is 
extended, with a scale length zvlba of approximately 
15 mas (compared to the angular scale of the orbital 
separation a of 0.1 mas), with an upper limit to the 
half-opening angle of ovlba < 2 ° where the viewing an - 
gle, 0LOS, is approximated to 40° (|Stirling et al.l[200lll . 
Combined with the other fiducial parameters of Cyg 
X-1, a robust limit can be derived from the observed 
stable compact jet by taking the bending angle to be 
V'oo = A^m^wind/^mjet < 7r/2, which gives 


L 


'jet 


1037 gpgg g-1 

r/(7 = 4/3) 


Sin 


^LOs \ 


-2/7 , , 

/aVLBA\ ^ 


/(t) 

T^wind 


Mwind 


sin 40° 

J V 2° 

Wlba 

N , n 2/7-1 

/150 

a 

/ J 

■yjet \ 

y 

.0.6c/ 

-1 


(ISOj^A-i 


(150)1/2 


2.6 X 10“° Mq yr 


r-i 


> 8.5 X 10 


-3 


( 20 ) 


where 7 = 4/3 in our standard model. For the case of 
7 = 5/3, the limit increases to3.9xl0“^ due to relatively 
shallow increase in jet thickness along the jet, which re¬ 
quires a larger initial opening angle oq and thus a wider 
initial jet to give the same observed ovlba ff, 2 °. 

This limit can be made more specific by the fact that 
the jet appears to be oriented in the same direction in 
separate VLBA radio observations, taken during different 
orbital phases (jStirling et al.ll2001^ . The data suggest a 
possible moderate bending on the jet of less then 10 ° on 
VLBA scales. Given the stability of the jet, we consider 
this an upper limit on the jet bending angle, i.e., ifoo ^ 
10 °, which translates to 


Ljet > 7.6 X 10^^ ergs s"i 
I Mwind 


2.6 X lO- 6 M 0 yr-i 


T'wind 


1.6 X 10® cms“i 


^ /aVLBA^ 

/ sindLCsV^^ 

/(b) 

/ V 2° / 

\ sin 40° / 

_/(7 = 4/3)_ 


/ 2:VLBA 

. V a 


) /150 


1-2/7 


(150) 


1-2/7 


(150)-V2 


( 21 ) 


where for 7 = 5/3, this limit increases to 3.5 x 
10 ®® ergss“®. 

This lower limit on the jet powe r is consistent wit h 
the range of jet im wers quoted in iGallo et al.l (|2005D . 
iRussell et al.l (|2007ll . and iSell et al.l ( 201511 . but derived 


from a completely independent, likely more robust argu¬ 
ment. 



Figure 9. Comparison of numerical results with analytic esti¬ 
mates. Each symbol represents the numerical result, and solid 
lines indicate the analytic jet trajectories. Dashed lines indicate 
the asymptotic towards which the jet is expected to converge ana¬ 
lytically. 



Figure 10. Fractional deviation of the jet trajectory between nu¬ 
merical result and analytic formula as a function of jet bending 
angle. 


4.5. Off-axis Jets 

When the progenitor of the compact object in an X-ray 
binary undergoes a supernova explosion in its last stage 
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of stellar evolution, it likely receives a substantial kick. 
Such a kick is capable of leaving the binary in an orbit 
where the s pin and the orbital angular rn omentum are 
misaligned (|Brandt fc Podsiadlowskilll995h . 

The spin will slowly align with the orbital axis through 
a combination of the Lense-Thirri ng effect and the inter¬ 
nal v is cosity of the a c cretio n disk ( Bardeen fc PettersonI 
I1975D . iMartin et al.l (|2008H showed that the alignment 
time scale is a few times smaller than the life time of the 
mass-transfer state, typically taHgn 10® —10®yr. While 
this implies that in most XRBs, spin and orbital angular 
momentum are likely aligned, it is plausible that they 
are mis-aligned in a sub-set of young XRBs, especially 
for HMXBs, given the short main sequence life time of 
the companion. 

In this situation, the jet-wind interaction will become 
phase dependent. Because the jet propagation and bend¬ 
ing time is short compared to the binary orbital period, 
the analysis presented above applies only to the phase of 
the binary orbit where the jet is perpendicular to the or¬ 
bital separation vector a. For any given orbit, there are 
two such node points. During other orbital phases, one 
side of the jet will approach the companion star (ie., the 
minimum distance between jet and star is smaller than 
the orbital separation), while the other side of the jet 
is receding; the approaching side will be more strongly 
bent. 

We performed a set of simulations to explore this sce¬ 
nario. We define the inclination angle 7r/2 — 0o between 
the initial direction of the jet z and the orbital separation 
vector a, such that 


9 q = 7r/2 — cos ^ 



( 22 ) 


i.e., 00 = 0 implies the jet is perpendicular to a, corre¬ 
sponding to the case discussed so far, while 0o = 7r/2 im¬ 
plies a jet maximally inclined, pointed at the star/away 
from it. In other words, in the limit of instantaneous 
reaction of the jet to orbital changes, the angle orbit 
between the jet and the orbital velocity vector does not 
affect jet bending (reflected also in the fact that most of 
our simulations neglect orbital motion entirely). Clearly, 
for any non-zero jet inclination relative to the orbital 
axis, the angle 0o will change as a function of orbital 
phase. 

Because jet bending reacts instantaneously to orbital 
changes, the jet will always propagate in a plane spawned 
by the initial jet direction and the orbital separation vec¬ 
tor, to order considered here. For circular binary orbits, 
the asymptotic bending angle therefore only depends on 
00 , and only implicitly depends on binary phase through 
00. 

We now extend our analytic formula to the off-axis 
case. The momentum flux of the jet will be the 
same, but the accumulated momentum flux of the stellar 
wind changes with inclination angle. The accumulated 
wind momentum per unit jet length is then (following 
eq. (IT4l)l: 


where 


f{z,l,do) = (cos0o) 


-2-bl/7 


/: 


dy 


1 


1 + 


(24) 

We can use the ratio of the accumulated wind momen¬ 
tum to jet momentum per unit length, eq. (HB, to derive 
the asymptotic bending angle tpoo 


where 


, Mwind^^windDet ~ 

V'oo = —7-r- - -/(7, do), 

47r a -Ljet,kin 

/(7,0o) = lim /(z,7,0o). 


(25) 

(26) 


It is obvious that if 0o is 0, this expression reduces to 
eq. (ITTI . 

For the approaching side of the jet, which is inclined to¬ 
wards the companion star in our simulation (Figure [TT]), 
/( 7 ,0o) is monotonically increasing, resulting in stronger 
bending of the jet. 

If 00 is larger than some value, the jet collides with the 
star. Given the binary parameters for Cygnus X-1, this 
limit angle is 62.18°. The impact region where the jet 
directly encounters the star is marked by the grey area 
in the Figure [Tljyill 

Figure dH shows 4>oc as a function of 0o, assuming hi 
is kept constant in the approaching side of the jet where 
it heads toward the companion star. From this analysis 
we can infer the jet bending angle for an initially off-axis 
jet. As one might expect, a more inclined jet will be 
more strongly bent. The right panel of Figure IT^ shows 
'4’ca with respect to the orbital plane. The figure shows 
a monotonic increase of jet bending angle with initial jet 
inclination. 

Figure [T51 shows the momentum ratio and asymptotic 
angle for the receding jet. For jets with higher inclination 
angles, the bending angle is reduced. 
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Figure 11. Variation of / as a function of 6o- The grey area 
shows the impact region where the approaching jet runs into the 
stellar surface. 
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Figure [14] shows density maps for several cases of the 
(23) misaligned simulations. All other parameters are the 
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top Side of the jets 
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namically curved jet. However, the jet bending we stud¬ 
ied in this work is likely not the cause of the phase lag, 
because the bending direction is parallel to the orbital 
separation vector, leading to only very small phase lags 
due to light-travel-time delays, while the observed phase 
lags would require bending in the direction against th e 
orbital velocity, as argued bv iSzostek fc Zdziarskil (1200711 . 

For higher jet inclination angles, {e.g. Sph- 
Wind_E36_60deg or SphWind_E36_75deg), both jet and 
counter-jet will impact the star once per orbit, leading 
to jet disruption and reformation; one might thus expect 
to observe episodic jet eruption from such a system. 


Figure 12. Case of approaching jet, initially pointing towards 
the companion star. Left panel: The ratio of the accumulated 
wind momentum flux to the jet momentum flux as a function of 
^0- Right panel: Asymptotic jet bending angle with respect to the 
orbital plane. The grey area shows the impact region where the 
approaching jet runs into the stellar surface. 


bottom side of the jets 




Figure 13. Case of receding jet, initially poi nting away from the 
companion star. Same panels as in Figure fl2l 

same as in our standard model, SphWind_E36. The 
inclination angles are 30°, 60°, 75°, respectively. The 
magenta solid line represents the analytic jet trajectory 
which is derived using the assumption of pressure bal¬ 
ance between the jet and the ambient medium. Eor all 
cases, the approaching jet shows some degree of disrup¬ 
tion and dynamical instability and the analytic approxi¬ 
mation breaks down for bending angles larger than about 
20°. On the other hand, the figures show that our ana¬ 
lytic approach is acceptable for the receding side of the 
jets. 

As one would expect from simple geometric consider¬ 
ations, while a jet launched perpendicular to the orbital 
separation a is bent symmetrically on either side, an in¬ 
clined jet shows asymmetric behavior between approach¬ 
ing jet and receding jet. Such a configuration would lead 
to increased jet bending in the approaching jet. Our 
limit on the jet power in Cyg X-1 was derived under the 
most conservative assumption that the jet is not inclined 
relative to the orbital axis. Because jet bending and dis¬ 
ruption increase for inclined jets, the possible inclination 
o f the jet will strengthen our limit on the jet power. 

iSzostek fc Zdziar^ ()2007l l argued that in Cygnus X- 
1 , the orbital modulation of radio emission occurs due to 
free-free absorption in the asymmetric wind as a func¬ 
tion of orbital phase, and the observed phase lag of the 
modulation with respect to the orbital phase could be at¬ 
tributed to the time delay for the emission from the dy- 


4.6. The case of Cygnus X-3 

We can apply a similar argument to the case of Cygnus 
X-3, an HMXB consisting of a compact object and a 
Wolf-Rayet companion. In this case, it is still not clear 
whether the compact object is a black hole or a neu¬ 
tron star, but the observational evidence of collimated 
j et-like structures is strong in multi-epoch radio maps 
(|Martf et al.ll200jT) . The radio jets with y-ray flares occur 
during bright soft X-ray states or during state transitions 
(ISzostek et al. [2008t Twani et al.ir2009fi . 

iDubus et al ( 2010 ) suggested that the jet of Cygnus 
X-3 is inclined relative to the orbital plane (20° < '0jet < 
80°) in order to obtain good fits to the gamma-ray modu¬ 
lation, with a significant offset between the jet footpoint 
and the site of the gamma-ray emission. 

This estimate does not take potential jet bending into 
account. In fact, as shown above, even if the orbital an¬ 
gular momentum and the jet axis are aligned, we should 
expect some bending to occur, and the asymptotic jet 
direction will change over the course of the binary orbit. 
As a result, the precession of the jet occurs on a time 
scale, Tprecess ~ Torbit in our model. We should point 
out that there is currently no evidence for (or against) 
jet precession in Cygnus X-3. The observed gamma-ray 
modulation can be explained by an inclined jet. An in¬ 
clined jet is likely bent more than a jet perpendicular 
to the orbital plane, so jet bending by the jet-wind in¬ 
teraction may have to be taken into account in detailed 
models of Cygnus X-3. 

Assuming the initial jet direction is perpendicular to 
the orbital plane, we can estimate the bending an¬ 
gle from eq. CZD from the binary and jet parameters 
Lifit = 10^®ergss“\ Mwind = lO'^Mgyr'^ a = 3 x 
cm, Uwin ri = 10 ^cms~^, and ujet = 0.5 c as used by 


det 

10 “ 


IDubus et ai.l (j20inH . 

Note that Cygnus X-3 is a considerably tighter system 
with higher inferred jet energy than Cygnus X-1, while 
the momentum flux from its neighboring Wolf-Rayet star 
is higher than that from OB companion in Cygnus X-1. 
In terms of the fiducial parameters given above, the jet 
bending angle can be calculated from eqs. (ED and (EZD 
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Figure 14. Density maps for off-axis jets, deviating from perpendicular direction to orbital plane: 30°, 60°, 75° from left to right panel. 
The solid magenta line indicates the analytical trajectory of the jet. The magenta circle indicate the location of the black hole. 
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indicating that for jet bending to be important, the ini¬ 
tial opening angle of the jet would have to be signifi¬ 
cantly larger than a few degrees, or the jet would have 
to be strongly misaligned with the orbital axis, with the 
jets very closely approaching the stellar surface for part 
of the orbit. 


4.7. Caveats 

Our simulations were carried out using an isotropic 
and homogeneous stellar wind model. This assumption 
allows us to study the dynamics of jet-wind interaction in 
a simple analytic framework. Several complications will 
affect this process in ways not included in this paper. 

4.7.1. X-ray Ionization 

Firstly, the X-ray flux from the accretion disk may be 
intense enough to ionize the circum-stellar gas, so the 
wind can be reduced or quenched, decreasing the jet 
bending angle. Line driving is likely inefficient in the 
illuminated portion of the wind if the ionization param¬ 
eter is above a critical threshold. Conversely, the X-ray 
illumination can itself produce thermal wind driving by 
heat input. The exact configuration of the wind is there¬ 
fore not entirely clear. 

However, from the fact that the black hole in Cygnus 
X-1 is accreting from the companion wind, the illumi¬ 
nated portion of the wind must be of similar density 
to that assumed in the shadowed region, which suggests 
that the estimates presented above will not be changed 
drastically by the effects of ionization of the wind. 

4.7.2. Clumping 


More importan tly, the stellar wind fr o m OB-type stars 
is like ly clumpy (lOskinova et al.l[2012tl . iPoutanen et al.l 
(|2008ll suggested that the wind in Cygnus X-1 is 
clumped, based on dips in the X-ray lightcurve. 

This situation was simulated by P12. Because global 
simulations that incorporate clumped winds with real¬ 
istic clump sizes and densities will be computationally 
impossible in the foreseeable future, given the large dy¬ 
namic range and resolution required to capture the mass 
in small clumps, our paper should be considered comple¬ 
mentary to the work presented in P12. 

While a detailed treatment of clumping in our simu¬ 
lations is beyond the scope of this paper, we will briefly 
discuss how significant clumping would affect the results 
of our simulations. For a fixed mass loss rate and wind 
velocity, the net effect of wind clumping will depend on 
the filling factor, the average size of the clumps, and the 
density contrast. Since secular bending is facilitated by 
the pressure gradient across the jet, we should expect 
bending in the presence of clumps to be reduced by the 
amount the momentum flux in the hot low-density back¬ 
ground component of the wind is reduced relative to an 
un-clumped wind. The reduction is set by the density 
ratio of the hot/low density background component phot 
to the density Pwind of an un-clumped wind of the same 
emission measure. 

Since phot is poorly constrained, we cannot easily quan¬ 
tify the reduction expected in the bending angle by the 
presence of clumps. Clearly, if phot is of the same or¬ 
der as the estimates of pwind used above, our results will 
be unaffected. However, if phot “C Pwind, the amount of 
secular bending of the jet will be reduced by the density 
ratio Phot/Pwind- 

In the case of 0-star winds in HMXBs, the presence 
of clumps in the illuminated side of the wind is even less 
well constrained, given the effects X-ray ionization may 
have on clump formation. 

For all these reasons, a general statement about how 
our results change in the presence of clumping is dif¬ 
ficult to make. However, considering the possible lim¬ 
its of clump sizes and density, we can estimate under 
which circumstances the jet would be able to propagate 
through the wind without disruption and derive similar 






























14 


constraints on the jet power for a given set of clump pa¬ 
rameters. 

If the density contrast between the clumps and the 
background wind is small, corrections should naturally 
be minor and the limit on the jet power from eq. 
will hold. If the density contrast is large and most of the 
mass is carried in clumps, the effect will depend on the 
average size of the clumps, relative to the width h of the 
jet: 

Large clumps: If the clumps are large compared to 
the cross section of the jet, jet-clump interaction will 
be disruptive, and the conclusions of P12 hold, i.e., 
the jet will be disrupted unless it is sufficiently power¬ 
ful to escape, in the case investigated by P12, Ljet > 
10 ^^ergss“^, similar to the conclusion we reach in this 
paper, with the main difference being that interaction 
with a clumpy wind int roduces significant sto chasticity. 
It has been suggested bv iZdziarski et aH (j20IIh that fluc¬ 
tuations in the radio light curve may be an expression of 
such interaction. 

Small clumps: If the clumps are much smaller than 
the cross section of the jet, such that the typical clump 
radius satisfies i?ciump ^ hjet, they will act like bullets 
passing through the jet, generating small bow shocks that 
slow down a small fraction of the jet. Small clumps will 
only globally disrupt the jet if (a) the covering fraction 
/clumps of clumps is larger than unity and (b) the total 
mass within clumps is sufficiently high to stop the jet 
fluid. Both conditions must be met for the jet to be 
disrupted by a wind whose mass flux is dominated by 
small clumps. 

Note that in this case, jet disruption will not be 
stochastic, since it requires the presence of many small 
clouds passing through the jet continuously, rather than 
a few large clouds. Since the jet in Cygnus X-1 is not 
disrupted continuously, at least one of the two conditions 
for jet disruption must be violated. We can now estimate 
the constraints the observed stability of the jet against 
wind disruption place on the jet and cloud parameters. 


must satisfy 


hi wind, clumped 


Pclump/vol ^^.^^ ^ vT^Mwind 

Pwind 


<Mwind 


I^clump 


< Mwind 


(30) 


(b) Lo-w mass density: A second necessary condition 
for clumps in the wind to significantly disrupt the 
jet is that the clump mass AMdump intercepted by 
the jet exceed the inertial mass of the jet plasma, 
Mjet. Conversely, a second sufficient but not nec¬ 
essary condition for the jet not to be disrupted is 

pOO 

^jet — / dzA^etp^et ,inertial 

^0 
poo 

I dzA^Qlfpcl■^^Ynp — AMdump (^1) 
which we can simplify as the condition 

Pjet ^ /voiPcluinp 


Assuming that the location of the re-collimation 
shock of the jet is set by the ram pressure of the hot, 
low density background gas of the clumped wind, 
we can use eqs. and ([8]) to write this condition 
as 

„,2 ^2 2 

n Pclump '^wind-' ‘jet.O ^^wind /QQ^ 

/vol- -2- = “2- 

^hot ^jet ^s,jet,0 

where phot "C Pwind “C Pdump is the density of the 
hot background wind in the clumped wind scenario. 
This condition can be written as a limit on the wind 
mass flux: 

2 

Mwind,clumped < Mwind < Mwind (34) 

Pwind Cg jet,0 


(a) Low covering fraction: The covering fraction 
/cover of clumps is given by 


where the last inequality reasonably assumes that 

Os,jet,0 ^ ^wind- 


dz 


f 


^/vol 


0 I^clump .^clump 


(28) 


In summary, we can distinguish the following three 
cases, one of which must apply to Cygnus X-1: 


where z is the path length along the jet, a is the 
orbital separation, and /voi is the clump volume fill¬ 
ing factor. If /cover ^ 1) the jet may be disrupted 
if the total mass within clumps is sufficiently high 
to stop the jet fluid (/dumps ^ 1 is a necessary con¬ 
dition for jet disruption). Conversely, a sufficient 
but not necessary condition for the jet not to be 
disrupted is that the covering fraction of clumps is 
small, /cover < 1, Or 


/vol « « 1 

a 


(29) 


The condition that the emission measure of the 
clumped wind is the same as that of the un¬ 
clumped wind (in the optically thin limit) gives a 
clump volume filling fraction of /voi = Pwind/Pciump 

^ fvo\ Pciump/Pwind — 1- ^Vith cq. (|29p , the con- 
dition on /cover implies that the wind mass loss rate 


1. The wind is not strongly clumped (that is, the mass 
flux in the diffuse/background wind is comparable 
to the mass flux used in this paper, even if clumps 
are present). In this case, our limit in eq. (I2I|) 
applies and the jet must be powerful to propagate 
through the wind. 

2. The wind is strongly clumped, with large clumps of 
size R> h] in this case, the analysis of P12 applies 
and the jet must be powerful to propagate through 
the wind. 

3. The wind is strongly clumped, with small clumps 
of size i? ^ /i; in this case, the mass loss rate must 
be orders of magnitude below the mass loss rate 
of an un-clumped wind, given the observed emis¬ 
sion measure and lack of jet disruption. In this 
case, we cannot use the observed lack of jet disrup¬ 
tion/bending to derive a limit on the jet power. 
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Thus, (a) either the wind mass loss rate must be orders of 
magnitude smaller than the nominal wind parameters for 
a uniform wind for the companion star, with important 
implications for wind formation in HMXBs, and/or (b) 
the jet power must be large [eq. d^ ]. similar to the 
conclusion reached from energy estimates of the large 
scale nebula. 

5. CONCLUSION 

We performed global hydrodynamic simulations to 
study jet-wind interaction in the wind of the companion 
star in HMXBs. The interaction results in the jet being 
bent with a characteristic bending angle that depends 
on the properties of the jet and the wind. Using the 
small bending angle approximation, we derived a simple 
analytic formula for the asymptotic bending angle. The 
formula is consistent with the numerical results, and it 
can be used to analyze observations of jets in HMXBs. 
We showed that the analysis is valid for bending angles 
smaller than about 20°. 

We applied the formula to two observed HMXBs, 
Cygnus X-1 and Cygnus X-3. We constrained the jet 
power in Cygnus X-1 to be Ljet ^ 10^® ergs s“^ from the 
lack of observed precession of the VLB A jet. This limit 
is consistent with previous estimates, but derived from 
completely independent arguments and likely more ro¬ 
bust. 

Given Cygnus X-3 parameters, we argued that the jet 
bending is likely not significant, unless the jet in Cygnus 
X-3 has a large opening angle a or is significantly less 
powerful than estimated in previous studies. 

The main caveat in the application of this model to ob¬ 
served HMXBs is the lack of knowledge about the prop¬ 
erties of wind clumping in massive stars. We discussed 
conditions under which clumping will affect the analysis 
of jet bending presented in this paper, and under which 
clumps themselves will lead to observable dynamic dis¬ 
ruption of the jet. We showed that, if wind clumping 
is dynamically important for the jet-wind interaction in 
Cygnus X-1, our limit on the jet power holds unless the 
wind mass loss rate is orders magnitude below the nomi¬ 
nal wind parameters derived for a uniform wind (in which 
case we would not be able to constrain the jet power from 
the lack of disruption or bending by the wind). 
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